arXiv: 1507.06648v2 [astro-ph.GA] 22 Sep 2015 


Astronomy & Astrophysics manuscript no. ms © ESO 2015 

September 23, 2015 


An extreme [O in] emitter at z = 3.2: a low metallicity Lyman 

continuum source 

S. de Barros 1 , E. Vanzella 1 , R. Amorm 2 , M. Castellano 2 , B. Siana 3 , A. Grazian 2 , H. Suh 4,8 ,1. Balestra 5 , C. Vignali 1,6 , 
A. Verhamme 7 , G. Zamorani 1 , M. Mignoli 1 , G. Hasinger 8 , A. Comastri 1 , L. Pentericci 2 , E. Perez-Montero 9 , 

A. Fontana 4 , M. Giavalisco 10 , and R. Gilli 1 

1 INAF-Osservatorio Astronomico di Bologna, via Ranzani 1,1-40127 Bologna, Italy, email: stephane.debarros@oabo.inaf.it 

2 INAF-Osservatorio Astronomico di Roma, via Frascati 33, 00040 Monteporzio, Italy 

3 Departement of Physics and Astronomy, University of California, Riverside, CA 92507, USA 

4 Harvard-Smithsonian Center for Astrophysics, Cambridge, MA 02138, USA 

5 INAF-Osservatorio Astronomico di Trieste, via G. B. Tiepolo 11, 34131, Trieste, Italy 

6 Dipartimento di Fisica e Astronomia, Universita degli Studi di Bologna, Viale Berti Pichat 6/2, 40127 Bologna, Italy 

7 Observatoire de Geneve, Universite de Geneve, Ch. des Maillettes 51, 1290 Versoix, Switzerland 

8 Institute for Astronomy, 2680 Woodlawn Drive, Honolulu, Hawaii 96822, USA 

9 Institute de Astroffsica de Andalucfa, CSIC, Apartado de correos 3004, 18080 Granada, Spain 

10 Astronomy Department, University of Massachusetts, Amherst, MA 01003, USA 

Received ; accepted 


ABSTRACT 

Aims. The cosmic reionization is an important process occurring in the early epochs of the Universe. However, because of observa¬ 
tional limitations due to the opacity of the intergalactic medium to Lyman continuum photons, the nature of ionizing sources is still 
not well constrained. While high-redshift star-forming galaxies are thought to be the main contributors to the ionizing background at 
z > 6, it is impossible to directly detect their ionizing emission. Therefore, looking at intermediate redshift analogues (z ~ 2-4) can 
provide useful hints about cosmic reionization. 

Methods. We investigate the physical properties of one of the best Lyman continuum emitter candidate at z = 3.212 found in the 
GOODS-S/CANDELS field with photometric coverage from U to MIPS 24pm band and VIMOS/VLT and MOSFIRE/Keck spec¬ 
troscopy. These observations allow us to derive physical properties such as stellar mass, star-formation rate, age of the stellar pop¬ 
ulation, dust attenuation, metallicity, and ionization parameter, and to determine how these parameters are related to the Lyman 
continuum emission. 

Results. Investigation of the UV spectrum confirms a direct spectroscopic detection of the Lyman continuum emission with S/N > 5. 
Non-zero Ly a flux at the systemic redshift and high Lyman-a escape fraction (f esc (Lycr) > 0.78) suggest a low H i column density. 
The weak C and Si low-ionization absorption lines are also consistent with a low covering fraction along the line of sight. The sub¬ 
solar abundances are consistent with a young and extreme starburst. The [O iii]/L14959, 5007+H/? equivalent width (EW) is one of 
the largest reported for a galaxy at z > 3 (EW([0 iii]/L14959, 5007 + H J3) ^ 1600A, rest-frame; 6700A observed-frame) and the NIR 
spectrum shows that this is mainly due to an extremely strong [O in] emission. The large observed [O iii]/[0 n] ratio (> 10) and 
high ionization parameter are consistent with prediction from photoionization models in case of a density-bounded nebula scenario. 
Furthermore, the EW([0 iii]/L14959 ,5007+H/7) is comparable to recent measurements reported at z ~ 7 - 9, in the reionization epoch. 
We also investigate the possibility of an AGN contribution to explain the ionizing emission but most of the AGN identification diag¬ 
nostics suggest that stellar emission dominates instead. 

Conclusions. This source is currently the first high-z example of a Lyman continuum emitter exhibiting indirect and direct evidences 
of a Lyman continuum leakage and having physical properties consistent with theoretical expectation from Lyman continuum emis¬ 
sion from a density-bounded nebula. A low H i column density, low covering fraction, compact star-formation activity, and a possible 
interaction/merging of two systems may contribute to the Lyman continuum photon leakage. 

Key words. Galaxies: high-redshift; Galaxies: evolution; Galaxies: ISM; Galaxies: starburst 


1. Introduction 

One of the most pressing, unanswered questions in cosmology 
and galaxy evolution is which are the sources responsible for the 
reionization of the intergalactic medium (IGM; [Robertson et al. 


et al. 


2015) and capable of keeping it ionized afterwards dFontanot 
2014; [[Giallongo et al.||2015| ). Until recently, it was diffi¬ 


cult to draw a consistent picture based on the different obser¬ 
vational constraints (e.g., reconciling the ionizing background 
and the galaxy UV luminosity density; |Bouwens et al.| 2015) . 
However, there is still a large number of unconstrained parame¬ 


ters at z > 6, the most important one probably being the Lyman 
continuum escape fraction (f esc (LyC); [Zackrisson et al.||2013| . 
Unfortunately, it is not possible to directly detect Lyman contin¬ 
uum emission at the epoch of reionization or even in the after- 
math (e.g., at any redshift z> 4.5), due to the very high cosmic 
opacity (e.g., Worseck et al. [2014 ). 

A number of surveys at 1 < z < 3.5 both from the ground 
and with the Hubble Space Telescope ( HST ), have looked for 
ionizing photons by means of imaging or spectroscopy and there 
have been some claims of detections (JSteidel et al.[2001[fShapley 
[eTaLl[20Q6l |Nestor et al.||20l3l |Mostardi et al.||2013| |2015p . 
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However, careful analysis of some claims combining different 
instruments (e.g., near-infrared spectroscopy and high spatial 
resolution probing LyC, [Siana et al.[ |2015| ) shows that the in¬ 
dividual detections must be considered with caution because of 
the difficulty to eliminate a possible foreground contamination 
( [Vanzella et al.|2010b[|Mostardi et al.|2015|). Several other sur- 


Similarly, it would be desirable to find Lyman continuum 
emitters at even higher redshift, within the first 2 Gyr after the 
Big-Bang (z > 3) and to look for properties and physical mech¬ 
anisms that allow ionizing photons to escape. |Amorin et al. 


veys reported on ly upper l imits (Siana et al.|2010[|Bridge et al 


2010 

Malkan et al. 2003; 

Vanzella et al. 2010b, 2012; Boutsia 

et al. 

2011; Grazian et al. 

2015). This could be due to the rar- 


ity of relatively bright ionizing sources, as a consequence of the 


( |2014a| ) report the properties of a z = 3.417 galaxy with possi¬ 
ble indirect signature of ionizing photon leaking such as a large 
[O iii]/[ 0 n] ratio, but only an upper limit on the Lyman con¬ 
tinuum emission. Very recently a class of candidate of strong 
H/3+[0 hi] emitters has been identified at much higher redshift 
(z > 6.5) by looking at the Spitzer /IRAC excess in channels 1 or 


combination of view-angle effects (Cen & Kimm 2015), stochas- 2 (depending on redshift; Schaerer & de Barros 2010; Ono et al. 

tic intergalactic opacity (Inoue & Iwata|2008 ( Inoue et al.{2014) 

2012) Labbe et al. 2013; Smit et al. 20141 2015; Roberts-Borsani 

and possibly intrinsically low escaping ionizing radiation on av- 

et al. 2015). The redshift of two of them have been spectroscopi- 


erage, in the luminosity regime explored so far (L > 0.1L*, e.g., 
[Vanzella et aH2010b| . 

At higher redshifts, the galaxy contribution to the cosmic 
ionization budget must be inferred from other properties that cor¬ 
relate with their ability to contribute to ionizing radiation, if such 
properties exist. This is an open line of research with ground and 
space-based facilities. Indirect signatures in the non-ionizing do¬ 
main have been investigated by performing photoionization and 
radiative transfer modeling. Low column density and covering 
fraction of neutral gas (that regulate the escape fraction of ioniz¬ 
ing radiation) are related to a series of possible indicators: 

1. the nebular conditions, (i.e. density-bounded or radiation 
bounded) t raced by peculiar line ratios like large [Q m]/[Q n] 
ratio (e.g., |Jaskot & Qey||2014 [Nakajima & Quchi|2014| ) 
or deficit of Baimer emission lines given the inferred star- 
formation_ rate (SFR) and the UV slope (F A oc T 3 ; 
|Zackrisson et al.|2013| ); 

2. strength of interstellar absorption lines as a trace o f neu¬ 
tral gas covering fraction (Jones et al.|2012||Heckman et al. 

Pont ; 

3. structure of the Ly a emission line, its width and position rel¬ 
ative to the systemic redshift (e.g., Verhamme et al.|20T5] ); 

4. morphology and spatial distribution of the star formation ac¬ 
tivity (e.g., nucleation) which seems to be another property 
that favors efficient feedback and eventuall y a transparent 
medium for ionizing radiation (e.g., [Heckman et al.pOTTj 
[Borthakur et al.|2014| ). 

While all these features are promising signatures, possibly us¬ 
able at the redshifts of the reionization, there is still a miss¬ 
ing direct link between them and the Lyman continuum emis¬ 
sion. In this respect there is an ongoing effort to investigate 
this possible link in a particular class of sources that seem 
to possess all the above properties. The idea is to look for 
“green pea” (GP) galaxies in the nearby Universe. GPs have 
been discovered in the Galaxy Zoo project (Cardam one et al.| 
2009). They are small compact galaxies in the SDSS survey at 
0.11 < z < 0.36, dominated by strong [O iii]/L14959, 5007 emis¬ 
sion (EW([0 iii]/U4959, 5007) > 100A, up to ~ lOOOA). This 
strong emission is thought to be related to a non zero ionizing 
photon escap e fraction (e.g, |Jaskot & Oey||2014{ |Nakajima &| 
Ouchi 2014). The typical physical properties of GPs are a low 
stellar mass (~ 10 8 5 M o < M* < 10 ia5 M o ), low metallicity 
(12 + log(0/H) = 8.05 ± 0.14, |Amorm et al |2010[( 2012|, high 


sSFR (sSFR=SFR/stellar mass, 1 to 10 Gyr \ Cardamone et al. 


2009), young age (few Myr), showing Wolf-Rayet signatures 


(e.g., bumps and He ii]/14686 emission line), and extremely large 
H a and [O iii]/15007 equivalent widths. Their strong Ly a emis¬ 
sion has been recently observed ( Henry et al.||2015| . However, 
currently, no direct LyC signatures have been identified in green 
peas and the work is still ongoing. 


cally confirmed at z = 7.73 and 8.67 with a combined equivalent 
width EW(H/3 + [O hi]) 720 A an d > 1 000 A (rest-frame), re¬ 
spectively ( [Oesch et al.|2015| [Roberts-Borsani et alJ2015; Zitrin 


|et al.|2015) . 

The nature of these objects is not well understood and further 
investigation is needed, though this will have to await the James 
Webb Space Telescope and the Extremely Large Telescopes. 

We have identified two Lyman continuum leakers in Vanzella 
et al. ( 2015| ), hereafter VI5: lonl and Ion2. These galaxies have 


been selected as Lyman continuum emitters through a photomet¬ 
ric selection which is based on the comparison between the ob¬ 
served photometric fluxes and colors probing the Lyman contin¬ 
uum emission and predictions from the combination of spectral 
synthesis models (e.g., Bruzual & Charlot|2003| and intergalac- 
tic medium (IGM) transmissions ( [Inoue et al.|20l4| . 

In this work we present the source with one of the largest 
[O ui]AA4959 ,5007 line equivalent width currently known 
(EW ~ 1500A, rest-frame) at z > 3 and with a plausible leak¬ 
age of ionizing radiation. Ultraviolet (VLT/VIMOS), optical and 
NIR (Keck/MOSFIRE) rest-frame spectroscopy are presented, 
as well as a detailed multi-frequency analysis with the aim to test 
the aforementioned signatures of linking Lyman continuum. We 
also use the extreme line ratios to investigate the nebular condi¬ 
tions. To this aim we also compare the observed Ion2 properties 
with theoretical predictions from photoionization models. 

Ion2 shares some of the properties of sources recently 
identified at much higher redshift (in particular the large 
EW(H/?+[0 hi]) and is therefore relevant to reionization, with 
the advantage to anticipate the rest-frame optical spectroscopic 
investigation of such extreme objects. 

The paper is structured as follows. The spectroscopic and 
photometric data are described in Sect. [2j In Sect. [3] we present 
the physical properties derived from the reanalyzed UV spec¬ 
trum and the newly acquired NIR MOSFIRE spectrum, and 
alongside we present the properties derived from the photom¬ 
etry. We also review previously reported properties (VI5) in the 
light of the new data. In Sect.|4j we discuss a possible AGN con¬ 
tribution to the observed Ion2 spectra. We summarize our main 
conclusion in Sect. \5\ 

We adopt a A-CDM cosmological model with Ho = 70 km 
s _1 Mpc -1 , £l m = 0.3 and £l\ = 0.7. We assume a Salpeter IMF 
(SalpeterI1955[). All magnitudes are expressed in the AB system 
(Qke & Gunn 1983). 


2. Data 

2.1. Spectroscopy 

The available low (LR, R=180) and medium (MR, R=580) res¬ 
olution VLT/VIMOS spectroscopy covering the spectral ranges 
3400-7000A and 4800-10000A have been obtained during 2005, 
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Fig. 1 . Postage stamps of Ion2 in . 6435 , V 60 6 , 1775 , Zsso, # 160 , IRAC 3.5/im, IRAC 4.5yum, IRAC 5.8yum, IRAC 8.0/im and MIPS 
24/tm. The sizes of the images are 4" x 4" for the HST bands and 20" x 20" for the Spitzer bands. 


with 4hr exposure time for both LR and MR spectra. The LR 
spectrum continuum sensitivity is S/N ~ 0.8 at 3700A with 4 
hours and m(AB)= 26.5 assuming seeing 0.8” and airmass 1.2. 
It is calculated over one pixels along dispersion and 2xFWHM 
along spatial direction. (Popesso et al. |2009| |Balestra et al.| 
2010 ). Slit widths of 1" were adopted and the median seeing dur¬ 
ing observations was 0.8". We have newly reduced both config¬ 
urations, in particular focusing on the Lyman continuum emis¬ 
sion (LR-spectrum) and the ultraviolet spectral properties in the 
non-ionizing domain (MR-spectrum). Reduction has been done 
as described in Bales tra et al.| (|2010|) as w ell as adopting the 
AB-BA method described in |Vanzella et al.| ( |201 l[|2014~b| ). Both 
methods produce consistent results. 

Recently, a near infrared spectrum of Ion2 has also been ob¬ 
tained with Keck/MOSFIRE during spring 2015 (PI: Gunther 
Hasinger). 32, 32 and 36 minutes integration time have been ob¬ 
tained in the /, H and K bands, respectively, with the AB-BA 
dithering pattern of 2.9". The adopted slit width was 0.7" pro¬ 
ducing a dispersion of 1.30A, 1.63A, 2.17A per pixel, in the J, 
H and K band, respectively; the pixels scale along the spatial 
direction was 0.18"/pix. Ion2 was inserted in a mask aimed to 
target also low-z objects. The J band is useful to monitor low-z 
emission lines. 


Reduction has been performed using the standard pipelin^] 
Briefly, two-dimensional sky-subtracted spectra, rectified and 
wavelength calibrated are produced, from which the one di¬ 
mensional variance weighted spectrum is extracted with a top- 
hat aperture. Wavele ngth s olution has been computed from sky 
emission lines (e.g., |Rousselot et al.|[2Q00| ) in the three bands 
reaching an RMS accuracy better than 1/20, 1/30 and 1/40 of 
the pixel size in /, H and K bands, respectively (typical r.m.s. 
< 0.05A). 


1 https://code.google.com/p/mosfire/ 


Flux calibration has been performed by observing the stan¬ 
dard star A0V HIP 26577. Given the slit width (0.7") and the 
median seeing conditions during observations (0.85"), particu¬ 
lar care must be devoted to slit losses. For this reason we check 
it by using a star included in the same science mask (WISE 
895), for which the continuum is very well detected in the entire 
MOSFIRE wavelength range, and whose /, H and K band to¬ 
tal magnitudes are known from CANDELS (Gro gin et al.|201l[ 
Koekem oer et al. ||20 11 ) and HUGS photometry (Fontana et al. 
2014| ) with typical S /N > 40. From the star WISE 895 and its 
near infrared photometry we derive a correction factor for slit 
losses and possible centering effects of about 50%. With this cor¬ 
rection the total line flux (H/3+[0 iii]/LI4959, 5007) derived from 
the MOSFIRE spectrum is fully consistent with the total line flux 
derived from the SED fit and K -band excess (see Section 13.3b. 
Ion2 is composed by two blobs separated by 0.2" (see Figure |l|f|) 
each one with a half-light radius Re - 0.3kpc, or less than 0.1" 
at z = 3.212. Therefore the slit losses for the entire system at 
the 0.7" MOSFIRE entrance slit are modulated primarily by the 
seeing during observations, which was on average 0.85". A cor¬ 
rection factor of 1.4 is consistent with the expected fraction of 
light falling outside the slit for a perfect centering of the target 
and Gaussian seeing with FWHM=0.85". The correction is also 
compatible with the discussion in |Steidel et al.[ ( |2014| ). Therefore 
a correction factor of 1.5 is reasonable in our configuration, ob¬ 
servational conditions and including a possible small centering 
uncertainty. 

Line fluxes, errors and limits from VIMOS and MOSFIRE 
spectra are reported in Table [l] 


2 Created with the Rainbow Navigator tool, https://rainbow. 
fis.ucm.es/Rainbow_navigator_public/ 
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Fig. 2. Two-dimensional LR VIMOS UV spectrum of Ion2 with different cuts to highlight the spectral features Lya and Ly/3 (bottom) 
and the Lyman continuum A < 912A, in the range 855-910A (top, vertical dotted blue lines). We show the moving average calculated 
within a rectangular aperture (855-910A x 1.25”, red-dotted rectangle) in the spatial direction divided by its r.m.s. on the left side. 
A signal is detected at A < 912A with S/N > 5. The inset shows the pixel distribution of the background after sky-subtraction in the 
region 855-910A (derived from the S/N spectrum). The distribution is symmetric ( skewness = -0.016) with the mean and median 
close to zero, +0.0057 and -0.014, respectively. No significant systematic effects are present. 


2.2. Photometry 


Deep U band imaging has been taken with the VLT/VIMOS in¬ 
strument providing a l-<x depth of m ag ~ 30 within 1 circular 
aperture. We refer the reader to |Nonino et ah (2009) f or a com- 


plete description of the data reduction (see also Vanzella et al. 

201 (B . 

The photometry from the U VIMOS band to the 
Spizter /IRAC 8.0/im band is taken from the CANDELS 


GOODS-S multi-wavelength catalog ( [Guo et al.|2013| ). The op¬ 
tical Hubble Space Telescope ( HST) images from the Advanced 


Camera for Surveys (ACS) combine t he da ta from |G iavalisco 
|et al.| ( |2004| ), |Beckwith et al.| ( [20061 ) and |Koekemoer et al 


(2011). The field was observed in the F435W, F606W, F775W, 
F814W and F850LP bands. Near-infrared WFC3/IR data com- 
bine s data from Grogin et al. ( 201 1) , [Koekemoe r et al. ( 201 1| ) 
and Windhorst et al.| (|2011|), with observations made in the 
F098M, F125W, and F160W bands. The V LT/HAWK-I K s band 
imaging are taken from the HUGS survey (Fontana et al. [|2014| ). 
The Spit zer /IRAC 1 and 2 data are taken from the SCANDELS 
survey ([Ashby et aL||2015| ), while the IRAC 3 and 4 are taken 
from [Ashby et al.[ (J2013J). Hereafter, we will refer to the fil¬ 
ters F435W, F606W F775W, F814W, F850LP, F098M, F125W, 
F160W as $435, V(506> hi 5 , ^814, Z850? ^098, J\i 5 > #160» respec¬ 


tively. 

We also check detection in Spitzer /MIPS bands and in the 
Herschel data. As shown in Figure [I] we are able to put an upper 
limit to the MIPS 24pm flux {S /N < 2, m^B = 22.27 ± 0.90, 
Santini et al. 2009[ ), while it is undetected in Herschel bands 
(< 1.2mJy, < 1.2mJy, < 2.4mJy inthe70yum lOOyurn, and 160/mi 
bands, respectively, 5 cr upper limits). 


3. Physical properties 

3.1. Spectroscopic detection and spatial distribution of the 
Lyman continuum emission 


Spectroscopic detection of stellar Lyman continuum emission 
at z ~ 3 has been previously claimed for few galaxies (e.g., 
|Steidel et al.||2001[ |Shapley et al. |2006| ) but this kind of detec¬ 
tion remains a difficult task , mainly because of foreground con¬ 
tamination ( [Vanzella et al.||2010bl |Siana et aL]|2015[ |Mostardi| 
et al. |2015] ), and a clear spectroscopic detection of Lyman con¬ 


tinuum emission at high-redshift is still missing. The only two 





UyiMOS ( <937A) _ACS/i775_HAWKI/Ks 


Fig. 3. Ground based VLT/VIMOS U-band observation of Ion2 
at a resolution of 0.2”/pixel (left panel). The contour of the sys¬ 
tem (derived from the ACS //775 band, center panel) is indicated 
and superimposed in blue to the VIMOS U-image and to the 
HAWK-I K s image (right panel). The cutouts are 2.6” vs. 2.6”. 


spectroscopi c detections reported in literature are from Steidel 


et al. ( 2001| ) in a stack of 29 LBGs and individually in a sub 


sequent spectroscopic follow-up, in 2 out of 14 LBGs observed 
by Shapley et al.| ( |2006| ), dubbed D3 and C49. | Vanzella et al.| 
( | 2010 b| ) performed dedicated simulations to deri ve the proba- 
bility of foreground contamination, specifically for Steidel et al.| 
( 2001]) and ( [Shapley et al. 2006) adopted observational setups. In 
both cases, the probability was high enough to rise the doubt on 
their reliability. Subsequently the two sources with possible LyC 
detection of Sha pley et al.| ( |2006| ) have not been confirmed as 
LyC sources. In particula r D3 was ruled out with much deeper 
narrow band imaging by Iwata et al. (2009) and |Nestor et al.| 
( |2011| ) and C49 was excluded from the list of LyC sources after 
the possible spectroscopic confirmation of a close foreground 
lower redshift source (Ne_stor et al. 2013| and finally confirmed 
to be a z=2.029 foreground source (Siana et al .|2015| ). Therefore, 
at present there are not direct LyC detections from star-forming 
galaxies, apart the source reported in the present work. 

The VIMOS LR spectrum of Ion2shows a clear signal de¬ 
tected blueward the Lyman limit. We performed a careful re¬ 
analysis of the Ion2 low-resolution UV spectrum by comput¬ 
ing the moving average of the flux at A < 910A and we find 
a signal with S/N > 5 (Figure [2]). While several systematics can 
affect the derivation of the LyC signal (e.g., background sub- 
Shapley et al.l|2006| ), we use different methods of sky 


straction 


subtraction (ABBA method and polynomial fitting with different 
orders) and they all provide consistent results. In particular the 
same moving average in the LyC region calculated over the S /N 
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Table 1. Emission line properties. 


Line 

Flux a 

Ly a 

15.4+ 1.0 e 

Ly a (blue) 

5.1 

Ly a (red) 

10.3 

N v/11240 

<0.16 

C iV/11550 

<0.18 

He 11 A 1640 s 

0.6 + 0.3 

O m]/ll666 g 

0.5 + 0.3 

C hi] /ll909 

1.3+ 0.3 

[O n]T3727 i 

1.5 + 1.5 

H/?s 

1.5+ 0.8 

[O m]T4959 

5.4+ 0.8 (7.4 + 1.1) 

[O hi] T5007 

22.1 + 0.8 


Dust corrected flux b 
15.4!f 0 ° 


< 0.4 

<0.4 

0 - 6 $ 

o.si 

l-3-o 

1.51J 

1 5 +1 
1 - J -o 

7.41’ 

99 i +1*0.8 
ZZ *-0.8 


FWHM C 
240 ± 60* 


400 + 90 h 


155 ± 70 
147 ± 30 


LW d z 

94 + 20 e 3.2183 1 + 0.003 

3.2126 ± 0.005 
3.2183 + 0.003 


18+* 3.2127 

< 50 3.2126 

112 + 60 3.2128 + 0.005 

365 + 70 3.2128 + 0.0008 

1103 + 60 3.2126 + 0.0006 


a Observed emission line flux and error corrected for slit los ses (50%) in unitsof 10 -17 erg s -1 cm -2 . In parenthesis, the [O m]/14959 flux 
adopting the theoretical ratio [O m]/15007/[0 m]/14959 = 3 ( Storey & Zeippen [2000j ). 
b The emission line dust correction method is described in Appendix | a|" 
c FWHM (instrumental corrected) and error in units of km s -1 . 
d Equivalent width and error in units of A. 
e Sum of the two components. 
f Main component. 

8 Marginal detections (Figure^: S/N (O m]T1666) ~ 1.7, S/N (He iU1640) ~ 2, and S/N(H/S) ~ 2 
h FWHM is the contribution of the two blended components, 1907A-1909A. 

1 The error is large because the line is on the top of a sky line (Figure [ 5 ]). 


two-dimensional spectrum we derived with ABBA method (see 
Vanzella et al. 2014a| produces stable results with no signifi¬ 


cant systematics and a signal (> 5<x) at the same spatial position 
of the target. This signal can be interpreted as a direct detec¬ 
tion of the Lyman continuum emission, but the presence of the 
two components in the HST images can cast some doubts about 
the origin of this detection. However, assuming that the faintest 
component keeps the same magnitude as the derived B 435 mag¬ 
nitude (27.25 + 0.24; VI5) at shorter wavelengths, the average 
S/N in the range 3600-3840A is ~ 0.5 per pixel (from ETC). 
Averaging over 20 pixels, as we did in Figure [2j we would ex¬ 
pect S/N ~ 2.2. In addition, the ZLband dropout signature of 
such a component (B-V = 0.62, VI5) between the #435 and V 606 
prevents a possible increased emission approaching the U- band, 
unless an emission feature is boosting the U -band magnitude. In 
such a case the only possible line would be [O ii]/13727 which 
would also imply a certain amount of star-formation activity, 
that in turns would be detectable through B aimer and/or Oxygen 
lines in the wide spectral range probed here. The most likely ex¬ 
planation is that the spectroscopic detection is due to a Lyman 
continuum emission emerging from the brightest Ion2 compo¬ 
nent. Furthermore, the ground-based VLT/VIMOS U -band spa¬ 
tial distribution shows that most of the U -band flux is emitted 
from the brightest component (Figure [5]). 


However, the U -band probes both ionizing and non-ionizing 
photons (A < 937A), so a fraction of the signal is not due only to 
Lyman continuum photons, and the ground-based observations 
are clearly limited in terms of resolution, as seen in Figure [3] 

The only way to clarify the exact position and the detailed 
spatial distribution of the Lyman continuum emission is to per¬ 
form dedicated HST observations. A proposal to observed Ion2 
with HST/F336W (17 orbits) has been recently approved (PI: 
Vanzella, cycle 23). Hopefully, emission line diagnostics can be 
used to characterize the gaseous and stellar content of Ion2 and 
provide some hints about a Lyman continuum leakage, as dis¬ 
cussed in the next Section. 


From the signal detected 
fesc(LyC) defined as: 

/ t it \intrinsic 
_ y^UV/^LyC) 

* esc ' rel ~ (F uv /F L ycY bserved 


at 880-910A, we derive the relative 

expose), (1) 


where (L uv / L Ly c) intrinsic is the intrinsic UV to ionizing luminos¬ 
ity density, (F uv / FL y c) observed is the observed UV to ionizing 
flux density, and exp (r Ly c) is the LyC attenuation due to the 
IGM. The signal measured in the LR spectum corresponds to 
~ 26.95 AB magnitude and misoo = 24.60 is the magnitude 
observed at 1500A rest-frame (using the V 606 band corrected 
for the Ly a flux). The average IGM transmission is 0.54 in the 
880 - 910A range ( Inoue & Iwata|2008| Inoue et al.|2014 ) and 
assuming Li 5 qq/L 9 qq = 3 ( [Nestor et al.|2013| ), we get: 

fescM = 0.64^}° ( 2 ) 


If we relax the assumption about the intrinsic £ 1500/^900 ratio, 
then f esc ,rei increases. Or conversely, assuming the maximum 
IGM transmission at z = 3.2 (0.74), Lisoo/^oo must be lower 
than 6.5 to keep f esc ,rei < 1- 


3.2. Physical properties of the interstellar medium 

As described in Section |2j we use VLT/VIMOS and 
Keck/MOSFIRE spectroscopy to obtain a wide Ion2 spectral 
coverage with 850A < A < 5700A (rest-frame). The emission 
line fluxes, full widths half maximum and equivalent widths are 
summarized in Table [T] Several lines are barely detected with 
signal-to-noise < 3 (e.g., H J3). However, thanks to the well con¬ 
strained redshift (see below) and because we know exactly the 
wavelength position of some features, we do not expect strong 
systematics. This is the case for the H/3 and [O ii]/13727 lines, 
whose signature are indeed at the expected positions. 

Individual line flux errors have been derived using the error 
spectrum. Additionally, we perform simple simulations to test 
the error associated to emission lines ratios in the MOSFIRE 
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Fig. 4. Two- and one-dimensional NIR spectra of Ion2 in the MOSFIRE H and /Cband. Two insets show the regions where H (3 
and the [O n]43727 doublet are expected accordingly with the redshift. The gray stripes in the one-dimensional spectra indicate the 
position of the sky lines. 


spectrum (Figure [4]). Figure [5] shows the positions where the 
[O hi] 45007 observed line has been added to the sky-subtracted 
spectrum (marked with numbers from 1 to 10) with different 
dimming factors (DIM=1, 3, 10, from top to bottom). This pro¬ 
vides a rough estimate of the typical error and S/N ratio of 
our measurements. For example, the typical error associated to 
the line ratio is [O iii]45007/[0 iii]44959 = 4.1 ± 0.7, which 
is consistent with the theoretical ratio (2.99) within 1.4-cr (if 
necessary, we fix the ratio to the theoretical one). Furthermore, 
the [O hi] 45007 line dimmed by a factor 10 (lowest panel in 
Figure [5} is similar to the observed Hj3 2D spectrum, indicating 
a likely [O iii]45007 to Hj3 ratio > 10. The marginal H (3 detection 
(with S /N ~ 2) allows us to derive aratio [O iii]/H J3 = 14.7+5.1, 
while for the [O n] line the non-detection translates to a 2 -cr 
lower limit of [O iii]45007/[0 n]43727 >10 (uncorrected for 
dust in both cases). The errors derived with our simulations are 
fully consistent with those derived from the error spectrum. 

The first property derived from emission lines is the spectro¬ 
scopic redshift. As already reported in V15, the C iii]41906.68 - 
1908.68 transition is clearly detected in the VIMOS MR spec¬ 
trum (with S/N = 8, Figure [6]). This feature shows a symmetric 
shape with a relatively large FWHM (= 400km s -1 ) with respect 
to other lines like Ly a and [O hi], suggesting the two compo¬ 


nents have similar intensities, even if they are blended and non 
resolved. The redshift we derive from C hi] is fully consistent 
with the redshift of Oxygen lines 4959-5007A identified in the 
MOSFIRE spectrum (see Table [T}. This provides a robust esti¬ 
mate of the systemic redshift, z = 3.2127 ± 0.0008. 


We derive ionization parameter, Oxygen and Carbon abun- 
dances using a modif ied version of the HII-CHI-mistry code 
( Perez-Monter o||2014 ), adapted to provide metallicity, C/O and 
ionization parameter in a Te-consistent framework, based on 
the comparison of the observed UV and optical nebular lines 
with a grid of cloudy photoionization models (Ferland et al. 


|2013 ). The lines used as input are C iv41550, C iii]41909, 
O iii]41666, H [J3, while the C iv41550/C iii]41909 (sensitive 
to the ionization parameter), C hi] 41909/0 iii]41666 (propor¬ 
tional to C/O), [O iii]45007/0 iii]41666 (sensitive to the elec¬ 
tron temperature) and (C iv41550+C m]41909)/Hy£> ratios are 
compared with model-based ratios in a ^-weighted minimiza¬ 
tion scheme. The derived abundances and ionization parameter 
are 12 + log(0/H) = 8.07 ± 0.44, log(C/0) = -0.80 ± 0.13, 
and log U = -2.25 ±0.81. We also derive these quantities us¬ 
ing all the lines except O iii]41666 (S /N ~ 1.7) and we obtain 
12 + log(0/H) = 7.79 ± 0.35 and log U = -2.32 + 0.11 (C/O 
can not be derived without O iii]41666). Both set of lines lead 
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Fig. 5. Observed and simulated two-dimensional MOSFIRE spectrum, with the observed [O iii]T5007 line replicated at different 
wavelength positions with different dimming factors (DIM). See text. 


to low metallicity and high ionization parameter (in the follow¬ 
ing, we consider the results obtain with all the line measure¬ 
ments and upper limits). Ion2 metallicity is similar to the typical 
green pea m etallicity (12 + log(0/H) = 8.05 ± 0.14, |Amorin| 
et al.|[2010] ). The metallicity and ionization parameter are also 
consistent with extreme emission-line galaxies up to z ~ 3.5 


(Amorin et al. 2014b|a| 2015| ). The metallicity of Ion2 is also 
consistent with the mean metallicity of star-forming galaxies se¬ 
lected through their extreme EW([Q hi]) ( [Maseda et al.||2014| 
Amorin et al. |2015 Ly et al. |2014[ ). Moreover, the Ion2 C/O 


line with a possible Lyman continuum leaking that could be ex¬ 
plained with a low neutral hydrogen column density. Also, the 
Ion2 extreme [O iii]/15007/[0 n] ratio, the metallicity, and the 
ionization parameter are consistent with cloudy models with a 
non zero Lyman continuum escape fraction ((Nakajima & Ouchi 
20 T 4 ] , Figure 11) 


ratio is strongly subsolar (C/O© = -0.26 [Asplund et aL||2 009) 
and consistent with the trend in C/O vs. O/H for local galax¬ 
ies (Garnett et al.[ [1999| ), which shows a nearly linear increase 
in C/O with Z >~ 0.2Z©. The low O/H and C/O abundances 
of Ion2 are similar to those found in other strong emission-line 
galax i es at z ~ 1-3 (e. g., Christensen et al.||2012[|Stark et al. 
|2014[ | James et ak]|2014| ) and suggest a significant contribution 
of O and C produced by massive stars, which is also consistent 
with very young and extreme starbursts. 

We compare the metallicit y and the ionization parameter 
with the results presented in Nakajima & Ouchi ( |2014| ) in 


Figure [7] Overall, Ion2 has a lower metallicity than all other 
galaxy populations presented in their work and one of the high¬ 
est ionization parameter. The higher ionization parameter is in 


The identification of Lyman continuum leakage from Green 
pea galaxies is a current line of research (e.g., jjaskot & Oey 
2014[ [Nakajima & Quchi||20T4} |Borthakur et al. |2()14[ |Yang 
et al.|2015| ) and the Green pea nature of Ion2 and its LyC leak¬ 
age represent the first concrete attempt to link these two prop¬ 
erties. Ion2 represents an extreme case of Green pea galaxy, 
being the highest redshift (z > 3) ultra-strong Oxygen emitter 
(with EW([0 m]/Ll4959,5007) ~ 1100A) currently known. The 
Lyman continuum leakage observed in Ion2 allow us to inves¬ 
tigate the relationship between the LyC leakage and physical 
and morphological properties. However, as shown in Stasiriska 
et al. ( |2015| ), the [O iii]/[ 0 n] ratio is also related to the spe¬ 
cific star formation rate and the metallicity: the [O iii]/[0 n] ra¬ 
tio increases with increasing sSFR and decreasing metallicity. 
For Ion2 the observed ratio is > 15 making Ion2 as an outlier 
in the Oxygen abundance vs. [O iii]/[ 0 n] ratio relation and the 
EW(H p) (i.e., sSFR) vs. [O iii]/[0 ii] ratio (Figure 2, [Stasinska 
et al.|2015[ >: Ion2 ratio is higher by ~ 0.6 dex for the observed 
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Fig. 6. ID and 2D UV VLT/VIMOS MR spectra. The transparent blue stripes indicate position of the emission sky lines. The 
transparent red stripes indicate the telluric molecular sky absorptions bands, B and A bands at ~ 6900A and ~ 7600A. In the same 
panel, three two-dimensional spectra are shown, the reduced, r.m.s. and signal to noise, from bottom to top, respectively. The inset 
on the left side shows the double peaked Lya line (from VI5). In the bottom-right a zoom in around He n/11640 and O m]Al666 
positions is shown. Their detection remain tentative. In the bottom-left, the ID spectrum is shown allowing to see the relative Ly a 
flux in comparison with other lines. 


EW(H J3) (Table[l| and the ratio is higher by ~ 1 dex compared to 
galaxies in the SDSS sample with similar metallicities. The Ion2 
[O iii]/[ 0 n] ratio is also higher tha n what is expected for th e 
derived stellar mass, SFR, and sSFR ( [Nakajima & Quchi|2014| ). 
Therefore, we conclude that the extreme [O m]/[0 n] ratio is due 
to unusual physical conditions (density-bounded nebula), which 
imply a low column density of neutral gas, and so favor leakage 
of ionizing photons ( [Nakajima & Ouchi[2014 ). 


3.3. Physical properties derived from the photometry 

While the emission lines highlight Ion2 ISM physical proper¬ 
ties, we also use the photometry to confirm findings such as 
the large [O hi] emission and derive other physical parame¬ 
ters. From Figure [8j the jump in the K-band is evident, K = 
22.97 ± 0.02 while the continuum level in the adjacent bands, 
/7 i6o and IRAC Channel 1 is - 24.15. The A m - 1.1 magni¬ 
tude in the HAWKI K-band corresponds to an additional flux of 
3.1 x 10 -16 erg/s/cm 2 , that corresponds to an observed equiva¬ 
lent width of 6700A (1600A rest-frame). The total flux derived 
for H/3+[0 m]/i/I4959,5007 from the MOSFIRE spectrum (see 
Tab. [TJ is fully compatible with the flux estimated from the K- 
band excess (after correcting for slit losses of the three lines, see 
above). The derived H/3+[ OIII] equivalent width is one of the 
largest ever observed at z > 2. The measured line ratios among 


the H/3+[0 m]/U4959,5007 suggest the [O m]/l5007 rest-frame 
equivalent width is ^ 1100A, with a similar value to the few ex- 
treme cases reported in literature at z ~ 1 - 2 ( Atek et al.|2011| 
Amorin et al.||2014b[ |Maseda et al.||2014|). The [Q m]/H ft r atio 


is also higher than the one reported in Holden et al.| (|2014) for 
LBGs at z ~ 3.5, while the ratio is higher in this sample than in 
the local Universe ( [Brinchmann et al.|2004| . Also at z > 6 large 
equivalent widths of H/3+[0 hi] have been reported, measured 
through photometric excess detected in the Spitzer/IRAC chan- 
nels (e.g, [Finkelstein et al.|2013[ [Smit et al.|2014[ [Oesch et al 


2015). While these cases can be spectroscopically investigated 


only with future facilities, like the James Webb Space Telescope 
(, JWST) and Extremely Large Telescope class, in our case the 
MOSFIRE near infrared spectroscopy has allowed us to access 
the details of the lines. 


We derive the SFR from the UV luminosity, H J3 and 
Lya lines (Kennicutt [1998J), neglecting dust attenuation (see 
Appendix |5J. We get SFR(UV) = 15.6 ± 1.5M 0 yr” 1 , 

SFR(H/5) = 40.0^ 5 5 6 M o yr 1 , and SFR(Lyof) = 14.3 ± 


1.0M© yr -1 . The three SFR indicators are consistent within 68% 
confidence levels, while SFR(H/3) is higher than the two oth¬ 
ers. The resonant nature of Lya photons prev ents us fro m de¬ 
riving precisely the SFR with this line (e.g., [Atek et al. 2014). 
On the other hand, the 1 1/3 line is barely detected (Table [I). The 
higher SFR derived from the H J3 line can be explained if the 
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Yang et al. 2015; Henry et al. 



Fig. 7. Comparison of the metallicity and ionization parameter 
(U = g/c) betw een the mean sample properties from Nakajima 
|& Quchi| ( |2014| ) and Ion2 with lcr uncertainties. Ion2 is com¬ 
pared to local galaxies from the SDSS survey, a galaxy sam¬ 
ple at z ~ 1, a Lyman break galaxy (LBG) sample at z ~ 2.5, 
a Lyman alpha emitter (LAE) sample at z ~ 2.5, local Lyman 
break analogu es (LBAs), LyC emitters, and Green Peas (GPs). 
See Nakajima & Ouchi ( 2014) and references therein for a com¬ 
plete description of these samples. 


age of the galaxy is significantly younger than lOOMyr because 
of the different equilibrium timescales between UV and nebu¬ 
lar emission lines. However, this is not the case according to 
the SED fitting, with an age for Ion2 ~ 200Myr. By comparing 
SFR(UV) and SFR(Lytr), we derive a Lyman-tr escape fraction 
f esc (\jy(i) > 0.78 (f esc (Lya) > 0.28 by using SFR(H/3) instead 
of SFR(UV)). This high Ly a escape fraction can be related to a 
blue [3 slope (e.g., Ha yes et al. | 2014|), but also to a low H i col- 


umn density (Ntj t < 10 iu cm 

|20B) . 

Ion2 Lya profile has been described in VI5: it exhibits a 
blue asymmetric tail. A low column density would be consis¬ 
tent with the bluer Ly a emission at the systemic redshift, and 
also with the high ionization parameter, and so with a leakage of 
ionizing photons. The star formation activity confined in a small 
physical volume can also favor the escape of ionizing radiation 
(e.g.,[Heckman et al.|2011 Borthakur et al.|2014| ). Furthermore, 
detailed analysis of the local blue compact dwarf galaxy NGC 
1705 shows that O stars are able to photoionize volumes far 
from their location ( [Annibali et al.||2015| }, which could create 
density-bounded conditions. Moreover, the spatial density of O 
stars in NGC 1705 increases toward the center of the galaxy, 
forming a sort of super star cluster efficient in ionizing its en¬ 
vironment. Ion2 could resemble a similar condition, being spa¬ 
tially resolved, but nucleated in the ultraviolet (1400A), in which 
a young super star cluster could be present in the core and signif¬ 
icantly contributing (or totally contributing if the AGN is absent) 
to the ionization field. Interestingly, because Ion2 has two com¬ 
ponents likely at a similar redshift (Section |L?| ), the leakage of 
ionizing photons can also be favored by the merging/interaction 
of these two components ( [Rauch et al.|2011| ). 

From SED fitting, we also derive a stellar mass 
log(M*/M 0 ) = 9.5 ± 0.2 and an age of the stellar population 
(since the onset of the star formation) log(Age/yr) = 8.6 ± 0.2. 


a (A) 

Fig. 8. Observed Ion2 SED (black and red dots). We show the 
best SED fit in blue and the integrated fluxes with blue crosses. 
We exclude photometric points (in red) affected by the IGM and 
strong emission lines (Lyor, [O hi]) from the SED fitting. The 
arrow indicates lcr upper limit. 


We notice that while the two Ion2 components are separated 
in the BViz bands, they are seen as a unique galaxy in the 
other bands, and the stellar mass is the total stellar mass of the 
two components. Given the very close separation among the 
two (0.2”), we decided to avoid any tentative decomposition 
and SED fitting of each component individually, especially in 
the Spitzer/IRAC bands where the pixel size is 1.22” and the 
decomposition is in practice not possible. Using SED fitting 
and anchoring the best fit to the ft slope derived from UV spec¬ 
trum fitting (Appendix [B|, we predict the expected observed 
L(IR)24/im assuming energy balance (i.e, we assume a balance 
between the energy absorbed in the UV/optical (we integrate 
over 912A to 3 jam) and re-radiated in the IR). The MIPS 24/tm 
observations probe the PAH feature at 6.2/mi (rest-frame) 
which can be used to derive the total IR luminosity (i.e., IR 
lumin osity integrated over 8-1000/im) and then t he bolometric 
SFR ( [Pope et aLl|2008[ |Rujopakam et al.||2012| ). Because the 
SFR derivation relies on several assumption (star formation 
history, age, metallicity, IMF, |Kennicutt||1998| ), we rely on the 
comparison between predicted L(IR)sed and derived L(IR) 24 yUm? 
rather than relying on a SFR comparison. From the MIPS 24yum 
upper limit, we derive log(L(IR)/L 0 ) 24 yUm = 9.92 +q^, while 
from SED fitting we derive log(L(IR)/L 0 ) SE D = 9.77^;™ 
The SED predicted and observed L(IR) are consistent but the 
uncertainties are too large to exclude an AGN contribution to 
the IR flux. We further discuss a possible AGN contribution in 
Section!?] 

Ion2 also shows an excess in the IRAC2 band, in compari¬ 
son with the IRAC1 channel. This excess can not be explained 
by emission line or other feature from both blobs (e.g., a contam¬ 
ination by a strong emission line as H a would imply z > 5.0 for 
the emitting galaxy). While we do not have a satisfactory expla¬ 
nation, we note that the rest-frame wavelength interval probed 
by IRAC2 (9500A-11900A) cannot be explored in detail be¬ 
fore the launch of JWST. Therefore that part of the spectrum 
is poorly known, especially for sources in the distant universe. 
In this work we do now comment further about this feature that, 
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however, do not affect our results. We summarize the physical 
properties of Ion2 in Table [2] 



-1000 0 +1000 
Velocity (km/s) 

Fig. 9. The velocity profiles of [O iii]/15007 and Ly a lines rela¬ 
tive to systemic. [O iii]/15007 line is shown in black, the Ly a first 
order in red, and the Ly a second order with the dashed blue line. 
The positions and width of the sky lines are marked with gray 
stripes. 


Table 2. Ion2 physical properties. 


log(M*/M 0 ) 

9.2 + 0.2 

log(age/yr) 

8.6+ 0.2 

SFR(M 0 yr‘) 

15.6 + 1.5 

E(B - V) stars 

<0.04 

12 + log(0/H) 

8.07 + 0.44 

log? 

8.23 + 0.81 

log(C/0) 

-0.80 + 0.13 


Interestingly, the Ly a profile exhibits a double peaked emis¬ 
sion, with the bluer component being half of the main red 
peak and having a redshift consistent with the systemic redshift 
(z = 3.212, Figure [9]). If the two Ly a components are due to 
the two Ion2 blobs, we would expect multiple-peaked [O hi] 
emissions ( [Kulas et al.|2012| ), which are not present in the Ion2 
MOSFIRE spectrum (see Section |T2| ). If the observed Ly a pho¬ 
tons arise only from the brightest Ion2 component, as discussed 
above i t can b e c onsi stent with a low neutral gas column density 
(e.g., |Verhamme et al.|2015|), and with the escape of Lyman con¬ 
tinuum photons ( Behrens et al.||2014| ). However, a small shift 
of the maximum of the profile compare to the systemic red¬ 
shift is expected in case of Lyman continuum emission (i.e., 
log(N H T ) < 18), with a maximum peak separation < 300km s 1 
( [Verhamme et al.|2015| . But the non detection of low-ionization 
absorption lines (e.g., [C n]/il334, [Si ii]T1260) is consistent 
with a low neutral hydrogen column density allowing the es¬ 
cape of ionizi ng photons and a low covering fraction along the 
line of sight ( [Heckman et al.|20lH| Jones et al.|2012||2013j ). We 
possibly detect O i + [Si n]/U303 with EW = -1.2 ± 0.7A and 
tentatively [Si ii]T1526 with EW = -1.0 ± 0.7A. However, the 
signal-to-noise ratio is low ( S /N ~ 2). Assuming a lcr fluctua¬ 
tion, we can put a lower limit on the EW of non detected lines as 
EW > -0.7A. This is in line with the results from Shapley et al. 
( 2003] ) at z ~ 3: galaxies with high Ly a equivalent width exhibit 
much weaker low-ionization absorption lines than galaxies with 
low EW(Lytr). 

Finally, Ion2 is not detected in the 6Ms X-ray Chandra im¬ 
age (L x < 2 - 3 x 10 42 erg s -1 ) and most of the high-ionization 
emission lines are not detected (e.g., N V/U240, C IV/U550). 
However, from the UV spectrum r eanal ysis, we report a tentative 
detection of He ii/!1640 (see Sect. |3.2| ): this line is is barely de¬ 
tected in the ABBA spectrum, but not in the Balestra et al. ( |2010| ) 
spectrum. Furthermore, VI5 also reported a possible variability 
of Ion2. We further discuss the possibility of an AGN contribu¬ 
tion in Section [4] 


3.4. Properties previously reported: summary 


Ion2 has been previously reported as a Lyman continuum emitter 
candidate in VI5. Here, we briefly summarize the Ion2 proper¬ 
ties reported in this latter reference. 

Ion2 is a compact source (300 ±70 parsecs) with two visible 
components. The presence of these two blobs can be misleading 
when trying to assess the detection of a possible Lyman con¬ 
tinuum emission, because multi ple components are often source 
of foreground contamination ([Vanzella et al. 2010b; Siana et al.| 
|2015[ |Mostardi et al.|[20l5] ). However, BViz magnitudes have 
been derived for each component after subtracting the close 
companion with a GALFIT fitting procedure (see V15). The two 
show the same break between the V^o6 and #435 bands and the 
absence of any spurious (low-z) emission line along the wide 


wavelength range probed here (from 3400A up to 24000A) sup¬ 
port that the faint component is a physical companion of the 
brightest one (i.e., at z = 3.2). Furthermore, the slit orientation 
in the MOSFIRE observation is also maximizing the separation 
in the spatial direction (see Figure [I]), though seeing FWHM is 
significantly larger than the separation. Therefore from photo¬ 
metric shape and spectroscopic arguments we conclude that the 
two components are at the same redshift. 


4. Is Ion2 an AGN? 


We now discuss the possibility that the main source of UV ioniz¬ 
ing photons is an AGN. First of all, the UV emission is resolved 
so it is most probably due to the stellar population, not an AGN. 
Therefore, an AGN could contribute to the SED but it would not 
dominate it at all wavelengths. 


Unfortunately, we can not use the BPT diagram ( [Baldwin 


|et al.[|1981| to discriminate between star-formation and AGN, 
because H a and [N 11 ] emission lines are out of the MOSFIRE 


wavelength range. Juneau et al. ( 2011| proposed an alternative 
to the BPT diagram, with the so-called MEx diagram based 
on the comparison between stellar mass and [O iii]T5007/H/3 
ratio. The probability for an object to be an AGN associated 
with the MEx diagram is not defined at the Ion2 position 
( Juneau et al.||2014j ), although the large [O m]/i5007 to Hfi ra¬ 
tio would classify Ion2 as an AGN. Other diagnostic diagrams 
have been proposed, like the [O iii]T5007/H/3 vs. [O n]T3727/H/3 
( Lamareille et al. 2004| ). Also in this diagram, Ion2 would be 
classified as an AGN, because of the large [O iii]T5007/H/3 ra¬ 
tio. However, for both diagnostics (MEx and [O iii]T5007/H/3 
vs. [O ii]/13727/H/3), there is no known AGN at a similar locus 
as Ion2. Alternatively, photoionization models predict a possi¬ 
ble increase of the [O iii]/[0 11] with decreasing metallicit y and 
increasing ionization parameter ( [Nakajima & Quchi|2014| ). 


10 
























































S. de Barros et al.: An extreme [O hi] emitter at z = 3.212 


As already stated, Ion2 remains undetected even in the 6Ms 
X-ray Chandra survey. According to the empirical relation be¬ 
tween the [O m]/i5007 luminosity and the X-ray luminosity (e.g., 
Pan ess a e t al.|2006] ), given the extreme [O iii]T5007 emission 
the expected X-ray luminosity would be L x ~ 10 45 erg s _1 if the 
[O iii]/15007 emission is associated to an AGN. Considering the 
scatter in the L([0 iii]T5007)-Lx relation, we can derive a lower 
limit of L x > 5 x 10 43 erg s _1 which would be easily detectable 
in the 6Ms Chandra data (with typical S /N > 10), independently 
on the source position in the X-ray image. However, extremely 
obscured AGN can have lower X-ray luminosity than expected 
from the L([0 iii]T5007)-Lx relation, this AGN possibly being 
Compton-thick AGN (e.g., Vignali et al.|2010| ). The infrared and 
X-ray luminosities derived from the MIPS detection and the 6Ms 
upper limit, respectively, can be c ompar ed with the properties of 
known obscured AGN (Lanzuisi et al. 2015[ ): Ion2 properties can 
be only consistent with a faint and highly obscured AGN (e.g., 
[Lutz et al.|2004| ). 

The large EW([0 iii]T5007) is unusually large for all types 
of AGN (Table [I] [Caccianiga & Severgnini|[2QTT ), while such 
large equivalent width has been alread y reported for some star¬ 
forming galaxies at lower redshift (e.g ., Atek et al.|2011[ van der 
|Wel et al.|2011t[Brammer et al.|2012 ). 

The UV spectrum is quite unexplored as a diagnostic for 
AGN/star-forming galaxies classification, but the observed nar¬ 
row lines (Table [I]) exclude the presence of broad-line AGN. 
We also compare Ion2 emission line ratios with typical UV 
emission line ratios in narrow line AGN. Nagao et al. ( 2006| 
present the UV spectrum analysis of narrow-line QSOs and 
narrow-line radio galaxies at 1.2 < z < 3.8 (see also De 
Breuck et ak |[2000| . The first evidence against a possible AGN 


contribution to Ion2 spectra is the lack of C ivT1550 detec¬ 
tion. This line is detected in all AGN types studied in |Nagao| 
et al. ( |2006| ). The typical line ratio C iii]/U909/C ivT1550 
is ~ 0.5 and C iv/11550/He n/11640 > 1, while we mea¬ 
sure > 10 and < 1 for Ion2 respectively. Narrow-line AGN 
UV spectra at z ~ 2-3 exhibit a C iii]/U909/C ivT1550 
ratio near 1 which is also much lower than the Ion2 ratio 
( [Hainline et al.||20TT] ). The C ivT1550/C m]/il909 ratio can 
be used to determine the nature of the ionizing source, with 
C ivTl550/C hi] Al 909 ~ 2 if the ionizing photons are produced 
by an AGN and C ivT1550/C m]/il909 ~ 0.001 if population 

2003| ). The clear 


I stars are the ionizing sources (Binette et al. 


C hi] Al 909 detection and the C ivT1550 non detection are con¬ 
sistent with a main ionizing source being stars rather than an 
AGN. 

The origin of He iiT1640 emission in galaxies has been mat¬ 
ter of investigation, possi bly produced by massive stars both in 
the lo cal Universe (e.g., Szecs i et"ak|2015[ Grafener & Vink 


2015 1 and at high-redshift (e .g., |Cassata et al.|2013||Sobral et al. 


2015 |Pallottini et al.||2015] ). He H/11640 detection is generally 


associated with AGN activity, but it has also been reported that 
He II/U640 associated with no C ivT1550 is a possible feature in 
star-forming galaxies at z ~ 3 ( [Cassata et al.|2013~] ). Thi s could 
be explaine d by a contribution from population III stars ( Sobral| 
et a l7||2015 ). However, as stated in the previous section, we re¬ 
mind that the He H/11640 detection is tentative. 

Summarizing, Ion2 is not detected in the 6Ms Chandra sur¬ 
vey and has no broad lines. The extreme [O iii]/H/3 ratio would 
lead to c onsider it as an AGN (Lamarei lle et al.||2004[ [Juneau] 
|et al.|2014| ), but this r atio can also be exp lained b y extreme ISM 
physical conditions ( Jaskot & Oey||MT3} |Nakajima & O uchi 


|2014| ). The la ck of X-ray detection can be explained with an ob 
scured AGN ( [Vignali et al.||2010| ), but the lack of C ivT1550 


detection in the UV spectrum is likely inconsistent with this hy¬ 
pothesis. Therefore, the main argument against an AGN contri¬ 
bution is the lack of the C iv/11550 emission line. If an AGN 
contributes to the Ion2 emission, it is likely a faint and obscured 


AGN with peculiar photoionization conditions (De Breuck et al. 

[ 2000 ] ). 


5. Conclusions 


In this paper, we present new observations with the 
Keck/MOSFIRE NIR spectrograph and a new analysis of the 
UV spectrum of a Lyman continuum emitter candidate. Ion2 
is an object at z = 3.212 composed of two distinct blobs, nu¬ 
cleated and resolved in the rest-frame UV, indicating emission 
from star-forming regions. Ion2 is also a candidate Lyman con¬ 
tinuum emitter with U band flux consistent with a leakage of 
ionizing photon (P(f eS c(LyC) = 0) < 10 -4 ), as reported in VI5. 
Contamination from foreground galaxies can be ruled out. 

Our main results can be summarized as follows: 


a new analysis of the UV spectrum shows a signal consistent 

with a direct detection of ionizing photons with S /N > 5 

and a relative Lyman continuum escape fraction f esc ,rei = 

0 64 +uo * 
u,u -0.10’ 

the Ly a emission at the systemic redshift, the high Ly a es¬ 
cape fraction, the non detection of low-ionization absorption 
lines are consistent with a low neutral hydrogen column den¬ 
sity, while velocity separation of the two Ly a peaks is in ten¬ 
sion with expectation (e.g.,[Verhamme et al. 2015| ; 
we find low metallicity (~ 1/6Z©), strongly subsolar C/O ra¬ 
tio and high ionization parameter (log U = -2.25) using a 
^-consistent method, in good agreement with previous re¬ 
sults at z ~ 2 - 3; 

Ion2 exhibit one of the largest [O iii]/[0 n] ratio observed at 
z > 3 and similar large ratios are predicted for galaxies with 


low metallicities and Lyman continuum leakage (Nakajima 
|& Quchi|2014] ); 

there is no clear evidence of AGN contribution to Ion2 emis¬ 
sion; 

the Ion2 physical properties (SFR, stellar mass, dust attenua¬ 
tion, strong emission lines) are similar to green peas, and the 
lon2 large EW(H/3+[0 iii]/U 4959, 5007) is similar with the 
ones derived for LBGs at z > 7. 


Ion2 physical properties, morphology, and very strong 
emission lines make it a possible analog of z ~ 8 LBGs. 
Furthermore, unlike local analogs (e.g., |Heckman et ak||2011| ), 
Ion2 probes Lyman continuum emitter properties only ~ 
1.5Gyr after the reionization epoch. From the very large 


EW([Q m]/L14959,5007+H/7) observed at high-redshif t (Oesch 


et al. 

2015 

Zitrin et al. 2015; Roberts-Borsani et al. 2015; Smit 

et al. 

2015 

), we can speculate that these galaxies have a leakage 


of ionizing photons. 


The direct detection of Lyman continuum emission leave lit¬ 
tle doubt about the fact that Ion2 emits ionizing flux. This galaxy 
certainly exhibits peculiar properties which however do not al¬ 
low to fully discriminate among different possible sources of 
ionizing radiation yet: stellar emission, faint and obscured AGN, 
or other ionizing sources. Only high resolution HST observations 
of the ionizing radiation may provide useful clues. In the near 
future, our approved proposal to observe Ion2 with //5T/F336W 
will hopefully shed new light on the nature of this source (PI: 
Vanzella). 
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Appendix A: Emission lines dust correction 

Emission lines are often dust corrected using the Baimer decre¬ 
ment (Hcr/H/5): the observed ratio between Baimer lines is 
compared with the expected ratio derived from theory (e.g., 
Ha/H/3 = 2.86; |Osterbrock]11989| ), and assuming an attenua¬ 
tion curve, the dust correction is derived at any wavelength (e.g., 


Kas hino et al.||2013[ Dominguez et al.||2013[ |Price et aL||2014| 

Steidel et al. 20 14j7 Unfortunately, regarding Ion2 , the H a line 

is outside the MOSFIRE wavelength range and we must rely 
on the dust attenuation derived from the UV p slope (fx °c AP\ 
Meurer et al. 1999) or from SED fitting (which basically also 
reli es on the fit of the U V /3 slope). 

[Calzetti et al.| ( |200()| derive a relation between the stellar and 
nebular color excess as, 


E(B - V)nebular “ 2.27 X E(B - V) stellar , 


(A.l) 


and this relation is often misunderstood as both color excesses 
being derived with the same attenuation curve. Actually the neb¬ 
ular color excess is derived from a line-of-sight attenuation curve 
(e.g., Cardel fiet al.|1989) , while the stellar continuum color ex¬ 
cess is derived with the Calzetti curve ( Calzetti et al. 2000]). 

This issue is challenged by contradictory result about the 
stellar color excess relative to the nebular color excess at z ~ 2: 
some studies found results consistent with the misinterpretation 
described previously (e.g., [Forster Schreiber et al.|2009|), while 
other found E(B - V) ne b u iar ~ E(B ~ V))steiia r (e.g., |Erb et al.' 
2006[ |Reddy et"aT]|2()l()[ Shivaei et al.||2015| ). A detailed com¬ 
parison of a large sample of z ~ 2 star-forming galaxies ob¬ 
served with MOSFIRE and with measured Baimer decrement 
for most of this sample ( [Reddy et al.||2015| ), shows that the re¬ 
lation between E(B - V) ne buiar and E(B - V) ste iiar is i n fact SFR 
(and sSFR) dependent (see also de Barros et al. 2015) . To derive 
the emission line attenuation, we use the relation provided by 
Reddy et aT7| ( [2015| ): 


E(B - V) nebular - E(B - V) stellar = -0.049 + 0.079 

x(log(sSFR S ED) + 10) 


(A.2) 


[Reddy et ah ( 2015) emphasize that Equation [A. 2| is affected by a 
large scatter (<x ~ 0.12). 


Appendix B: UV p slope and SED fitting procedure 

We want to derive the UV p slope as accurately as possible be¬ 
cause this quantity is known to be related to the dust attenuation, 
while also dependent on the age of the stellar population and 
the star formation history (Leitherer & Heckman 1995} Meurer] 
|et al.|1999) . To derive the UV (3 slope, we perform a multi-band 


fitting, using the bands between 1500A and 2500A (Castellano 
|et al.|2012) . However, the flux in the zsso band is lower than what 


we expect from the observed fluxes in the band just blueward and 
redward zsso (hu and To 9 s)- We compare theoretical expectation 
from the MOSDEF attenuation curve and the MW curve ( Seaton 
1979) with Ion2 observed colors in Figure [Bn] The MW curve 


exhibits the 2175A bump which could explain the flux decrease 
in the 50 band, while uncertainties are too large to discrimi¬ 
nate between a curve with or without the 2175A bump. However, 
Scoville et al. (2015) report that their data are consistent with the 


0.6 


0.4 


0.2 


o 

LO 

00 


10 

CM 


0.0 


- 0.2 


-0.4 


- 0 . 


': 8 . 


Increasing 



MOSDEF attenuation curve (Reddy et al. 2015) 
MW attenuation curve (Seaton 1979) 


75 


-0.5 


-0.25 0.0 

F775-F125 


0.25 


0.5 


Fig. B.l. Range of F775-F125 and F125-F850 col ors pred icted 
with BC03 templates and t he MOSDEF ( [Reddy et al.|2015| ) and 
MW attenuation curves (Seaton 1979). Ion2 colors are shown 
with the red square. 
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Fig.B.2. ID UV spectrum from V15 (black). ACS photometry 
(blue squares) is shown along with the corresponding transmis¬ 
sion curves (V 606 , ins, hu, and Zsso, blue dashed lines). We also 
show the medium band photometry from the MUSYC survey 
(green stars, [Cardamone et al.||2010) . The red dots and associ¬ 
ated error bars s how the flux derived through the wavelength 
windows used in [Calzetti et al. ( 1994) to derive the UV p slope. 
The cyan line denotes the best fit over the Calzetti windows 
and the violet and yellow lines encompass the 68% uncertainty 
((3 = -2.75, -3.40 and -2.20, respectively). 


presence of the 2175A bump, while for Reddy et al. (2015) this 
presence is only marginal. 

We derive the /3 slope from the UV spectrum, using the wave¬ 
length windows described in [Calzetti et al.[ ( l994[ ), but exclud¬ 
ing the 10th window because being outside the UV spectrum 
wavelength coverage. The median flux within each window is 
calculated and used to fit a linear relation in the observed spec¬ 
trum, assuming the flux density per u nit wa velength Fx ~ AP, in 
units of erg s" 1 cm -2 A -1 ( Hathi et al. 2015). The resulting slope 
is p = -2.75^-^. It is worth noting that the UV slope from the 


spectrum is not influenced by the possible 2175A bump. The UV 
slope derived from the spectrum is consistent within ler with 
the p slope derived from multi-band fitting of the photometry 
(- 2.20 ± 0 . 20 ). 
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Reg arding the SEP fitting procedure, we use BC03 tem¬ 
plates ( [Bruzual & Charlot|2003| and we anchor the result to the 
UV [3 slope by not taking into account solutions inconsistent with 
the p slope derived from the UV spectrum, within l<x. We also 
fix the metallicity to Z = 0.2ZO, which is consistent with the 
metallicity derived in Section |3.2| We test different star forma¬ 
tion histories (declining, rising, constant) and attenuationcurves 
(Calzetti, MOSDEF, MW), assuming a Salpeter IMF ( |Salpeter| 
1955). The results show that the physical parameters such as the 


stellar mass and the age are not strongly affected by the assump¬ 
tions. Therefore, in this work, we assume a constant star forma¬ 
tion history and the MOSDEF attenuation curve. While the p 
slope is consistent with no or low dust attenuation, we use the 
relations derived in |Reddy et al. ( 2015| ) to translate the p slope 
in color excess. Thisjeads to E(B - V) ste ii a r < 0.04, while the 
Meurer relation (Meur eret al.|19 99) leads to no dust attenuation 
because of the higher intrinsic p slope in this latter relation. 
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